Multifractal behavior of linear polymers in disordered media 
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The scaling behavior of linear polymers in disordered media modelled by self-avoiding random 
walks (SAWs) on the backbone of two- and three-dimensional percolation clusters at their critical 
concentrations p c is studied. All possible SAW configurations of N steps on a single backbone 
configuration are enumerated exactly. We find that the moments of order q of the total number of 
5^4 Ws obtained by averaging over many backbone configurations display multifractal behavior, i.e. 
different moments are dominated by different subsets of the backbone. This leads to generalized 
coordination numbers [i q and enhancement exponents j q , which depend on q. Our numerical results 
suggest that the relation (ii = p c pL between the first moment /xi and its regular lattice counterpart 
[i is valid. 
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I. INTRODUCTION 



The question of how linear polymers behave in a disordered medium has attracted much attention in recent years. 
The problem is not only interesting from a theoretical point of view, but may also be relevant for understanding 
transport properties of polymeric chains in porous media, such as in enhanced oil recovery, gel electrophoresis, gel 
permeation chromatography, etc. In this context, it is useful to learn about the static or conformational 

properties of linear chains, modelled by self-avoiding walks (SAWs), in the presence of quenched disorder, e.g. how the 
surrounding structural disorder influences their spatial configuration. As a quite general model of a random medium, 
percolation may be considered as the paradigm for a broad class of disordered systems and has therefore been 
mostly used so far. 

We are interested in how the statistical behavior of SAWs on percolation clusters at criticality (p — p c ) differs 
from their behavior on regular lattices. While the values of the exponents for SAWs on regular lattices are well 
established pp|-^2]], there is no complete agreement about their values on percolation clusters at p c |l3|,[l4j]. Here, we 
study (i) the so-called effective coordination number of the cluster, where contradicting results have been reported 
using different numerical techniques. Next, we consider (ii) the enhancement exponent 7 and (in) the exponents v r 
and vi characterising the end-to-end distance of SAWs in the r- and l-space metrics. Finally, we determine (iv) the 
values of the critical exponents describing the corresponding structural distribution functions. 

We concentrate on SAWs on percolation clusters at p c in two and three dimensions. In the literature, two distinct 
methods have been used for evaluating 571 Ws: Exact enumeration (EE) and Monte Carlo (MC) simulation. In the 
EE technique, all SAW configurations on a given cluster are taken into account, but only relatively short chains 
can be evaluated. In a MC simulation, longer chains can be studied, but inherently the ensemble of configurations 
remains incomplete. Here, we use the EE technique in combination with an appropriate finite-size scaling procedure 
to determine the relevant exponents. Since 'infinitely' long chains can only exist on the backbone of the cluster, 
where dangling ends are absent on all length scales, we study the SAWs directly on the backbone. This enables us to 
generate longer chains on a given cluster and to average over a larger set of different cluster configurations. 

Specifically, we enumerate all possible SA W configurations of TV steps for a single backbone and study different 
moments of the total number of SA Ws and of their end-to-end distance by averaging over many different backbone 
configurations. Our analysis shows that the critical exponents v r and vg do not depend on the order q of the moments, 
while the enhancement exponents and the effective coordination numbers do depend on q, leading to multifractal 
behavior. In particular we find that the first moment of the effective coordination number fix satisfies fix = p c [i, 
where fj, is the effective coordination number of the underlying regular lattice, resolving previous controversies. The 
mean structural distribution functions for the end-to-end distance after A steps, both in Euclidean and topological 
space, are obtained numerically, supporting the expected scaling forms Jl5|Jl6| ]. 

The paper is organized as follows: In Section [n[ we briefly review the main relevant properties of SAWs on regular 
lattices to illustrate the different numerical procedures employed in this work. In Section III, we present results for 
the total number and the mean end-to-end distance of SA Ws on the backbone of the incipient percolation cluster. 
The corresponding distribution functions of the end-to-end distance and their scaling behavior, in Euclidean and 
topological space, are also discussed. Finally, in Section IV, we summarize our main results. 



II. SAWs ON REGULAR LATTICES REVIEWED 



In this section we illustrate the different numerical techniques we use in the following by briefly reviewing the main 
results for SAWs on regular lattices. The main idea is to show that our finite-size scaling employed in the later sections 
enables us to obtain quite accurate estimates for the critical exponents based on EE results for relatively short chains. 
Here, we consider the case d = 2, which is particularly suitable since many results are known exactly. 



A. Total number of SAW configurations Cjv 



The total number Cat of SA W configurations of A steps behaves as |fl]] 

Cjy = A fi N TV- 1 (1) 

where fi is the effective coordination number of the lattice, 7 is the universal enhancement exponent, and A is a 
constant. To determine /i, 7, and A, we choose to study the behavior of the quantity 

In Cat In A , . „, In A 
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as a function of N. Fig. |l] shows for the square lattice, that the values for \i and 7 obtained by fitting the EE data 
using Eq. (0) agree well with the accepted values reported in literature (see Table |) . 



B. Mean end-to-end distance and structural distribution function 

The root mean square end-to-end distance of SA Ws of N steps, r(N) = [^(N)] 1 / 2 , averaged over all possible SAW 
configurations behaves as 

r(N) oc N V¥ (3) 
with the universal exponent vp — 3/4 in d — 2 as suggested by Flory M . In Fig. |[ we show values for r(N) versus 



N obtained by EE technique 
as shown in the inset of Fig. 



|2| an 



The asymptotic value for vp (see also Table |) is obtained using successive slopes 
and is in excellent agreement with the theoretical prediction. 
A more detailed information about the spatial structure of SA Ws is given by the distribution function P(r,N), 
where P(r, N) dr is the probability that after N steps the end-to-end distance of a chain is between r and r + d?\ 
This quantity obeys the scaling form [ fulfill 

P(r 7 N) oc —f(r/N VF ) (4) 
r 

and is normalized according to dr P(r, N) = 1. The analytic form of the scaling function f(x) is known asymp- 
totically: 

f(x) OC I X x g 2+ d exp (_ cx 5j ^ I J \ , (5) 

where g\ = (7 — 1)/^f # 2 = $ (d Wf — 1/2] — [7 — 1]) (22) and S = 1/(1 — vp) (!(]]. Values for these exponents are 
summarized in Table Q. We have verified these predictions by enumerating all SA W configurations for N = 23 and 24 
and calculating the corresponding distributions P(r,N), from which we have extracted the different exponents (see 
Fig. ||) . We show that a more accurate determination of the exponent <? 2 compared to a simple fit using Eq. (^J) can 
be obtained by employing a specific numerical procedure described in Appendix |A] (see inset of Fig. ||). The obtained 
values are in agreement with the theoretical predictions (see Table ||) . 



III. SAWs ON THE BACKBONE OF THE INCIPIENT PERCOLATION CLUSTER 

Next, we consider SA Ws on the incipient percolation cluster by generating all 5^4 W configurations directly on the 
backbone of the cluster. We obtain the backbone of a given cluster grown by the Leath algorithm |2^,|24| by randomly 
choosing one of the sites of the last grown cluster shell (e.g. site A in Fig. and determining the backbone between 
site A and the seed of the cluster (site S in Fig. ||) by the burning procedure described in p5p^ ]. The SAWs start 
at the seed S of the cluster. To avoid finite size effects, the chemical distance between both endpoints S and A of 
the backbone is chosen at least 20 times larger than the chemical length of the SAWs. The large ratio between both 
chemical lengths is needed, since close to the endpoint A, the backbone has a quasi-linear structure, which would 
falsify the results for the 5^4 Ws. The straightforward idea to use all sites on the last grown shell as endpoints for the 
backbone does not help, but introduces boundary effects in the opposite direction, since in this case the backbone 
coincides with the cluster near the end points; cf. j26]. 

We analyse the results for SAWs on the incipient percolation cluster by applying analogous numerical procedures on 
the data as described above for SAWs on regular lattices. In contrast to the case of regular lattices, on a percolation 
cluster two different metrics can be defined, the Euclidean metric and the topological metric. On average, the chemical 
distance I between two backbone sites separated by the Euclidean distance r increases with r as 

I oc r dml " (6) 

where d min = 1.1306 ± 0.0003 in d = 2 || and d miD = 1.374 ± 0.004 in d = 3 ||]. Thus Eq. (§) yields the scaling 
relation between the two metrics, which will be used in what follows. Numerically it is found, that data obtained 
in ^-space show less fluctuations (cf. e.g. |fL5[ ). Therefore more accurate estimates for many characteristic quantities 
(such as critical exponents) in r-space can be determined by studying the corresponding quantity in ^-space and 
transforming it to r-space. For example, the fractal dimension of the backbone in £-space is df — 1.45 ± 0.01 in d = 2 
and df = 1.36 ± 0.02 in d = 3. Using Eq. (§), this leads to the values df = dfd min = 1.64 ± 0.02 and df = 1.87 ± 0.03 
in r-space, respectively pq]. 
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A. Total number of SAW configurations: Multifractality 

Due to the disordered structure of the clusters, the total number Cn.b of SAW configurations that are generated 
on a single backbone, with the seed S of the cluster as starting point, fluctuates strongly among different backbone 

configurations. To characterize these fluctuations, we study the moments (c%&J- A similar study on percolation 

clusters at criticality has been performed for 'ideal' chains, i.e. chains which can intersect themselves. This model 
leads to a non-trivial dependence on q [|TJ. 
In generalizing Eq. (Q), we make the ansatz 

(c% >B ) 1/q = A q ^N^- 1 , (7) 

where fiq are the generalized effective coordination numbers of the backbone and j q the generalized enhancement 
exponents. Results for different values of q are shown in Figs, ^(a) and |^(b) for the square and the simple cubic 
lattice, respectively, employing the numerical procedure described in Section [0]A The values for [i q and "f q are 
displayed in Fig. |5| for d — 2, clearly revealing a dependence on q, reminiscent of a multifractal behavior. For large 
negative values of q, backbone configurations with a small number of SA W configurations Cn.b are singled out in the 
averaging procedure. We find that [i q — > 1 and j q — ► 1 for q — > -co, pointing to rare configurations of backbones 
with an almost linear shape. On the contrary, for large values of q the averaging procedure emphasizes backbone 
configurations with a large number of SAW configurations Cn,b- Since these backbones are the most compact ones, 
fi q and 7 9 are strongly enlarged. Fig. || seems to suggest that the structure of the most compact backbones differs 
distinctively from the structure of a regular square lattice, as lim^oo /i q « 1.9, which is well below the value for /i on 
the regular square lattice, and lim^oo j q w 1.7 is well above the value for 7 on the regular square lattice. 

These results resolve earlier controversies regarding the values for both [i and 7 for percolation obtained from MC 
simulations and by EE techniques. For the square lattice, for example, the values /i porc (EE) = 1.53 ± 0.05 jj2| and 
7porc(EE) = 1.33 ± 0.02 |53| have been obtained from exact enumerations calculations, while from MC simulations 
the values Mpcrc(MC) = 1.459 ±0.003 and 7pcrc (MC) = 1.31 ±0.03 @ were determined. We find fjti = 1.565 ±0.005, 
71 = 1.34±0.05 and fi = 1.456±0.005, 70 = 1.26±0.05, corresponding to the EE and MC results, respectively. This 
can be understood by noting that EE calculations yield by definition the whole ensemble (the so-called 'annealed' 
average), corresponding to the case q = 1, i.e. the normal arithmetic average. In contrast, MC simulations intrinsicly 
sample only a small subset of all possible configurations, omitting rare configurations, yielding 'typical' subsets of the 
ensemble (the so-called 'quenched' average). This quenched average is usually described by a logarithmic average, i.e. 

(Cjv,B) typ = CX P (InCw.B), and is equivalent to the limit q — > of Eq. ([?]), i.e. lim 9 ^ (^n = CX P ( m @n,b)- 

Indeed, our results are in excellent agreement, in both d = 2 and d = 3, with the relation 

Hi = PcM (8) 

where \i is the effective coordination number of the underlying regular lattice, p c — 0.5927460 for the square lattice |3?| 
and p c = 0.311605 for the simple cubic lattice p6| . This relation, which was originally suggested in the form 
Mpcrc = PcM lHI' could not be confirmed earlier on because of the different values obtained for /i pcr c. 

Because of the possible existence of rare events playing a dominant role in the average procedure, we have per- 
formed a detailed analysis of our numerical data to confirm that we have considered a sufficiently large set of cluster 
configurations (cf. Appendix |B|) . 



B. Mean end-to-end distances and structural distribution functions 



Next, we study the scaling behavior of the distribution functions for the end-to-end distance, (Pb((,N)) and 
(Pb(t, N)), averaged over many backbone configurations, where Pb(£> N) d£ is the probability that after N steps the 
chemical end-to-end distance of a chain on a single backbone is between £ and £ + d£, and ffe(r, N) dr is the analogous 
quantity in r-space. These distribution functions are expected to obey scaling forms similar to the one valid on regular 
lattices, Eq. (^), with the corresponding scaling exponents jlE|. The mean chemical end-to-end distance (1(N)) and 

1/2 \ 



the root mean square Euclidean end-to-end distance (r(N)) 



'-(N) 



scale with TV as 



(£(N)) oc N u 



(9) 
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and 



(r(N)) oc N" r , (10) 

respectively. The first average is performed over all SAW configurations on a single backbone, the second average 
is carried out over many backbone configurations. Following Eq. (Q), the exponents vi and v r are related to each 
other b y v r = z^/d m i n . The numerical results for and v r obtained by the successive slopes technique discussed in 
Section II B for regular lattices are reported in Table |n[ As an example Fig. |^ shows the determination of vi in d = 3. 



Accordingly, the scaling variable in chemical space is £/N v '- and the mean structural distribution function, averaged 
over many backbone configurations, has the form 



with the scaling function 



(PB(£,N))K~f(e/N>") , (ii) 



x gi+df exp (- Cd lX S ?) , X > 1 

Equivalently, in r-space, the scaling variable is r/N Ur , and one has 



(P B (r,N)) X -f(r/N^) (13) 



r 



with 



M*)<*\ l^i' _ / . ^ ! " ! • (14) 



x 



ff J+df cxp ^_ CdrX 5 r ^ ^ 



Both distribution functions are normalized according to J °° d£ (P B (£,N)) — 1 and / °° dr (Pb{t\N)) = 1- 

The numerical results for the distribution functions in d — 2 and d = 3 are shown in Figs, p] and ^, respectively, 
in both I- and r-space. The values for the exponents vi and v r = fi/d m i n reported in Table □ are used in the 
scaling variables. For the determination of the exponents g{ , g^, <?i, and according to Eqs. (|l2 ) and flli|) , we use 
the previously reported values of the fractal dimensions and c?P p6[ . The exponents and g[ can be estimated 
directly from the slope of fi and f r in the double logarithmic plots. Since g[ and g\ are related by g[ = g\ d m i n |fL5| , a 
more precise estimate for g\ can be derived from the estimate for g\. The determination of g\ and g r 2 is more difficult, 
since both exponents occur in the non-dominant part and are masked by the exponential. Therefore it requires the use 
of the slightly more involved numerical procedure discussed in Appendix [a| (see the insets of Figs. || and || for d = 2 
and d = 3, respectively). The numerical results we obtain for g\, g^, g[ and g^ are reported in Table |l[ Regarding 
the exponential factors, our results for the exponents Se and 8 r are consistent, within the present accuracy, with the 
expressions S( — 1/(1 — v i) an d 5 r = 1/(1 — v r ), respectively. 



As discussed in Section II B , for regular lattices the exponents g\ , v-$ and 7 are related by the des Cloizeaux relation 
gi = (7 — l)/vp. Therefore, it is legitimate to ask if a similar 'generalized des Cloizeaux' relation holds also for 
SA Ws in percolation. Since the enhancement exponent j q depends on q, it is necessary to find out whether the 
exponents vt and g\ as well as v r and g\ depend on q. To this end we generalize the averages (£{N)) and (f(N)) to 

" q (N)y /Q oc N^ q) and (r q (N)) 1/q oc N"™. Since this is equivalent of studying the quantities [/ t q P B (£, N) d£] 1/q 

and [J ri P B {r, N) dr] 1/q ', respectively, and (P B (£,N)) and (P B {r,N)) scale with t/N vt and r/N v % also ^Z 9 (iV)^ 1/9 

and (r q (N)) 1 l q must scale with £/N Ve and r/N Ur , respectively. Therefore vf^ — vt and vf^ = v r for all q, which we 
confirmed numerically. Furthermore we have numerically verified that the exponents gf and g\ (as well as g%, g\ 1 &i 
and 5 r ) are independent of q. Regarding the 'generalized des Cloizeaux' relation, our numerical results suggest that 
in d = 2 the relations g± = (7^=1 — 1)/^ and g\ — (7^=1 — l)/fr hold, where the, to some extent arbitrary, choice of 
7 g= i is motivated by the fact that q = 1 describes the annealed case of the whole SA W ensemble. However in d = 3 
these relations are not satisfied by the present numerical results. 
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IV. CONCLUDING REMARKS 



We have studied structural properties of SAWs on the backbone of the incipient percolation cluster in both d = 2 
and d = 3, applying exact enumeration techniques. Our results suggest that 5^4 Ws display multifractal behavior, 
caused by the underlying multiplicative process yielding an infinite hierarchy of generalized coordination numbers 

fi q and enhancement exponents j q describing the moments (cjj of the total number of SAWs of length N. The 

present results resolve previous inconsistencies regarding the suggested relation /x pGlc = p c /J>, where p c is the percolation 
threshold of a specific regular lattice, and \x and /K per c are the corresponding effective coordination numbers of SAWs 
for the ordered case and on the incipient percolation cluster, respectively. We find that this relation is accurately 
obeyed on the square and simple cubic lattice by identifying /i pcic = H\. 
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APPENDIX A: IMPROVED PROCEDURE FOR THE DETERMINATION OF g 2 



The procedure used for extracting the exponents <?2, g\ > an< ^ 92> describing the scaling form of the structural 
distribution functions, is an improved version of the procedure by Wittkop et al. |37j (cf. Q ) and is illustrated here 
for the case of regular lattices. The distribution function Eq. (0) can be written as 



O R 

P(r,N) = —f(r/N^) 
r 

with Q = 2n in d = 2 and Q = An in d = 3 and the scaling function f(x) defined as 

u gi+d , X < z 



/(*) = 



x 92+d exp (-bx s ) , x > z 



(Al) 



(A2) 



where 8 = 1/(1 — Vp). The actual value of the crossover z is determined from the numerical results. The constants B 
and b can be obtained from the normalization condition 



and from the second moment 



P(r,N)dr 



r 2 P(r, N) dr = r 2 (N) = iV 2 ^ 



Upon integration of Eqs. (A3) and (A4), one gets the exact relations 



B 



1 
9. 



1 



-r 



-,bz' 



?gi+d 



gi + d 



where T(m, z) is the incomplete Gamma function, and 

1 



ilB 



Sb ( -92+d+2)/8 



r 



9 1± d+2 - 



z- 



gi+d+2 



gi + d + 2 



(A3) 



(A4) 



(A5) 



(A6) 



Thus, by plotting the distribution function in the case x > z as y = b {92+d ^ 5 (Q.B)- 1 r P(r, N) exp (b l / 5 r/N v *) 

versus b 1 ^ r/N UF in a double logarithmic plot, t he e xponent gi can be read off from the relation y ~ x 92+d and 
adjusted until the above relations Eqs. (|A5|) and (A6) are satisfied. This method yields much more accurate results 



G 



than by directly fitting the distribution function itself. The accuracy of the result can be assessed by plotting 

y = -Inlb^+V^iClB^r P(r, N) (b 1 / 6 r/N UF )~ i92+d) ] = b (r/N u *) s versus b l ' 5 r/N Vv in a double logarithmic 

plot, from which the exponent 8 can be determined and compared with the expected value 5 = 1/(1 — up). The 
procedure can be extended straightforwardly to study the distribution functions {Pb{^,N)) and (Pg(r,N)) of SAWs 
on the backbone of critical percolation clusters. 



APPENDIX B: GENERALIZED AVERAGING PROCEDURE 



To obtain an estimate of whether the ensemble B of backbone configurations considered is sufficiently large in order 
to get convergent results, we analyse the data by a generalized averaging procedure as follows: The total ensemble B 
containing n to t backbone configurations is divided into subsets Bj containing n e g configurations each. The generalized 
average is then defined as 

Mii^Efe)?) " . (ED 

The obtained results (Cjv.b)^ depend sensitively on the different backbone configurations and display strong fluc- 
tuations, indicating that the system is not self-averaging. In order to smooth out these fluctuations, a second average 
is performed. This second step is a logarithmic average over the n to t / n c g subsets |39| 



C 



N,B 



(q, n cS ) = exp (in (C N , B )^ t ) = A w ^ {{ N^' 1 . (B2) 



In Eq. (B2), the limiting case n e g = 1 corresponds to the limit q — ► 0, while the usual average (cf. Eq. (]?])) is 
recovered when n c g = ntot- The results for the coordination numbers (J>q,n e a & n d enhancement exponents 7 g ,n eff 
are shown in Figs. |T^(a) (for q = 1) and |l^(b) (for q — 2). A dependence of these two values on n c g indicates 
that the given ensemble is too small to obtain the asymptotic values. If, on the contrary, the ensemble of backbone 
configurations is sufficiently large, then /i 9l n off and 7 q ,n off no longer depend on n c g. For q = 1, this seems to be the 
case when n e g > 10 3 , and for q = 2 when n c ff > 10 4 . 
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literatu&/32 a 


pt fl-.-Hl fcsgiflfe 




2.6385 ±0.0001 b 


2.641 ± 0.005 


u v 


3/4 c 


0.745 ± 0.005 


9i 


ll/24 d 


0.4 ±0.1 




5/8 c 


0.61 ±0.05 


8 


4* 


4.5 ±0.5 



TABLE I. Structural parameters for SAWs on regular lattices in d = 2. Results of t he pre sent simulations obtained on the 
square lattice, compared with the accepted values from the literature: a Ref. @, b Ref. c Ref. §], d Ref. c Ref. Q, 

and f Ref. @. 



7i 


1.34 ±i(H)2 


1.29 ±Z(M)8 


7o 


1.26 ±0.05 


1.19 ±0.05 


Mi 


1.565 ± 0.005 


1.462 ± 0.005 


Mo 


1.456 ± 0.005 


1.317 ±0.005 


v t 


0.89 ±0.01 


0.910 ± 0.005 


v r (directly from data) 


0.778 ±0.015 


0.66 ±0.01 


V r = Vll dmin 


0.787 ±0.010 


0.662 ± 0.006 


g{ 


0.45 ±0.10 


0.66 ±0.15 


9l = 9l dmin 


0.51 ±0.11 


0.91 ±0.20 




1.6 ±0.16 


1.95 ±0.17 


92 


1.26 ±0.18 


2.96 ±0.18 




9.5 ±0.5 


12 ±0.5 


S r 


4.85 ±0.20 


3.1 ±0.2 



TABLE II. Structural parameters for 5^4 Ws on the backbone of percolation clusters at criticality in d = 2 and d = 3, on the 
square and simple cubic lattice, respectively. The values for v r obtained directly from the numerical data are in agreement with 
the more precise values obtained from the relation v r — v^jd^^. The values for g[ = g 1 d m i n are also in good agreement with the 
corresponding values obtained directly from the data. The numerical values for the exponents g\ and have been determined 
using the procedure described in Appendix |X| Note that there is no simple relation between g\ and g^, i.e. g^ 7^ g| d m in- The 
values of 8e and 8 r are consistent, within the present accuracy, with the expressions Se — 1/(1 — and 5 r — 1/(1 — u r ). 







FIG. 1. The total number Cjv of SAW configurations on the square lattice plotted as (lnCiv)/JV versus N, from the presently 
available exact enumeration results for Cm, N < 51 |l7j]. The continuous line corresponds to a numerical fit obtained in the 
range 10 < N < 51 using Eq. (|), with (j, = 2.641, 7 = 1.3 and A = 1.35. 



FIG. 2. The mean end-to-end distance f(N) versus N for SAVFs on the square lattice. The continuous line is drawn as a 
guide and its slope has the theoretical value up = 3/4. In the inset, the successive slopes i/p = d lnr(iV)/d In N are plotted 
versus 1/N. A linear extrapolation of the points to the limit 1/N — > yields our estimate z^f = 0.745 ± 0.005, consistent with 
the value 3/4. 



FIG. 3. The structural distribution function of SAWs, rP(r,N) versus r/N VF , with vy = 3/4, for N = 23 (diamonds) and 
N = 24 (circles) on the square lattice. The dashed line in the range r/jV" F < 1 has a slope gi + d = 2.4, and the one for 
r/N" F > 1 is a fit with Eq. (|), for x > 1, yielding g 2 +d = 2.9±0.4, S = 4.5±0.8, and c = 0.7±0.1. In the inset, we show the 



function rP(r,N) = b {s2+d)/s (VLB)- 1 r P(r, N) 



exp 



= b (r/N" F ) s versus b 1 ^ 6 V /N" F , following the procedure 



(b 1/5 r/N v * 

described in Appendix |aJ allowing a more precise determination of g 2 . For our estimate of the crossover value z = 0.4 the 
continuous line has a slope g 2 + d — 2.61 ± 0.05, in agreement with the theoretical value (see Table F 



FIG. 4. A percolation cluster on the square lattice (full squares) and its corresponding backbone between the seed S and a 
site A randomly chosen on the last grown shell. 



FIG. 5. Generalized moments (C N B ^ of the total number Cn,b of 5^4 W configurations on the backbone of critical percolation 



clusters, plotted as (1/N) In 



<^, B ) 1/9 



versus N. (a) On the square lattice, for q = 2, 1, 0.5, 0, —0.5, —1 and —2 (from top to 



bottom); (b) on the simple cubic lattice, for q — 1 (top) and q = (bottom). Averages over 10 a backbone configurations each 
are performed. The continuous lines are the best fits based on Eq. (pj) , yielding the values for fi q and j q for q = and 1 given in 
Table [ll]. Some representative values for 7 9 , in addition to those reported in Table [n], are: 7-2 = 1.15 ±0.05, 7_i = 1.23 ±0.05, 
and 72 = 1.36 ± 0.05 in d — 2. Values of A q are found to fluctuate in the range of 1.0 to 1.3 in both d = 2 and d = 3. 



FIG. 6. The effective coordination numbers fi q and enhancement exponents 7, versus q for — 10 < q < 10 in d = 2 obtained 
from Fig. ^(a) . Except for j q for q > 2 the errorbars are smaller than the symbolsizes. The values for it and 7 on regular square 
lattice are marked by arrows, clearly showing that limg-^oo 7, is larger than 7 on regular square lattice. The inset shows fi q 
versus q for —2 < q < 2 in d = 2, in good agreement with the theoretical result fi q = /Uo(l + go"o/2) (continuous line), expected 
for \q\ -* JTeJ, with ^0 = 1-456 and a = 0.45. 



FIG. 7. The mean topological end-to-end distance (J(N)\ versus TV for 5^4 Ws on the backbone of critical percolation clusters 
in d = 3 averaged over 5 • 10 4 backbone configurations. In the inset, the successive slopes vi = dln(i(N)^ /dlniV are plotted 
versus 1/N. A linear extrapolation of the points to the limit 1/N —* yields our estimate ui — 0.910 ± 0.005. 



FIG. 8. Scaling plots of the distribution functions on the backbone in d = 2, for N = 39 and 40, averaged over 5 • 10 3 config- 
urations, (a) I (Pb(£,N)) versus £/N" e : The dashed line has the slope 1.90 and corresponds to the ansatz Eq. (|l^), for i<1; 
the continuous line is a fit with the ansatz Eq. (p^), for i>l, yielding gi = 1.4±0.4, Si = 9.5±0.5, and c 2 ,i = 0.09±0.01. The 

inset shows £(P B (t,N)) = b { f 2+d * )/Se (JIB^ 1 i (P B (£, N)) exp (h^^l/N 1 '^ = b t p"') f ' versus b^ &l t/N v >- , with our 

estimate of the crossover value zt = 0.21, according to the procedure described in Appendix [a|, yielding the more precise estimate 
g 2 +df — 3.05±0.15 (continuous line), (b) r (Pb (r, N)) versus r/N Vr : The dashed line has the slope 2.15 and corresponds to the 
ansatz Eq. (p^), for 1 < 1; the continuous line is a fit with the ansatz Eq. (p^[), for x S> 1, yielding = 1.46±0.4, S r = 4.9±0.3, 

and c 2 , r = 0.79 ±0.10. The inset shows r (P B (r, N)) = b { f 2+d ^ )/Sr (VlBr)" 1 r (P B (r,N))exp (bl /Sr r /N^A^ = b r (r/N^) Sr 

versus b Y J &T r/N" T with our estimate of the crossover value z r = 0.25, according to the procedure described in Appendix |X| 
yielding the more precise estimate g 2 ± dj 3 = 2.9 ± 0.15 (continuous line). 
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FIG. 9. Scaling plots of the distribution functions on the backbone in d = 3, for N = 39 and 40, averaged over 5 ■ 10 3 config- 
urations, (a) £ (Pb(£,N)) versus i/N" 1 : The dashed line has the slope 2.02 and corresponds to the ansatz Eq. (p^), for 
the continuous line is a fit with the ansatz Eq. (|l^) , for x > 1, yielding gi = 1.3±0.6, S e = 12.0±0.5, and c 3 , £ = 0.06±0.01. The 

inset shows £{P B (£,N)) = b ( / 2+d * )/Se (SlBe)' 1 £ (P B (£,N))exp (b\ ,h £/N v ^j St = b t (£/N" e ) Se versus b\ /5l £/N" 1 , with our 

estimate of the crossover value zt = 0.4, according to the procedure described in Appendix^, yielding the more precise estimate 
gi + df =3.31±0.15 (continuous line), (b) r (P B {r,N)) versus r/N" r : The dashed line has the slope 2.78 and corresponds to 
the ansatz Eq. (fl4|), i < 1; the continuous line is a fit with the ansatz Eq. (|l4|), i>1, yielding = 2.3 ± 0.6, S r = 3.5 ± 0.5, 

and c 3 , r = 0.88 ±0.10. The inset shows r (P B (r, N)) = bi 9 * +d ™ )/Sr (SIBr)- 1 r (P B (r, TV)) exp 



versus b^ Sr r /N Vr with our estimate of the crossover value z r 
(continuous line). 



(bl /Sr r/N v 

0.5, yielding the more precise estimate g\ + d 



b r (r/N"*) Si 
= 4.83 ±0.15 



FIG. 10. The effective coordination numbers p g ,n ef( (circles) and the enhancement exponents 7 9 ,n cff (squares) of SA Ws on 
the backbone in d = 2 for: (a) q = 1 and (b) q = 2 versus the effective ensemble size n e ff- The values are obtained by 
a least-square-fit of (\n[CN, B (q,n e ff)]) /N = (\n A q: n et! ) / N + In fJ.q,n eS + [(7<j,n off — 1) \nN]/N versus N, shown as fJ,q t n et c and 
7<j."ctf ~ 1 versus n e ff. 
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Ordemann et al, Fig. 5 
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